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Abstract 

The impact of heavy mediators on neutrino oscihations is typically described by non-standard 



■ four-fermion interactions (NSIs) or non-unitarity (NU). We focus on leptonic dimension-six effective 

. operators which do not produce charged lepton flavor violation. These operators lead to particular 

r~| \ correlations among neutrino production, propagation, and detection non-standard effects. We 



point out that these NSIs and NU phenomenologically lead, in fact, to very similar effects for a 
neutrino factory, for completely different fundamental reasons. We discuss how the parameters 
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I and probabilities are related in this case, and compare the sensitivities. We demonstrate that the 
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■ matter NSIs is broken for hadronic neutrino production, while NU is a fundamental effect present 

at any experiment. 



NSIs and NU can, in principle, be distinguished for large enough effects at the example of non- 
standard effects in the ^-r-sector, which basically corresponds to differentiating between scalars 
and fermions as heavy mediators as leading order effect. However, we find that a near detector at 
superbeams could provide very synergistic information, since the correlation between source and 
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I. INTRODUCTION 



During the past decade, experimental studies of neutrino oscillations have provided us 
with compelling evidence that neutrinos are massive particles and lepton flavors mix. Since 
an important new window is opened for searching new physics beyond the Standard Model 
(SM) of particle physics, it is interesting to discuss the impact of potential non-standard 
effects^ on neutrino oscillations. In this study, we focus our attention on non-standard effects 
from heavy mediators, which are integrated out at the scales of the neutrino oscillation 
experiments. 

It is convenient to parameterize the impact of the heavy flelds, present in high-energy 
theory, by adding a tower of effective operators O'^ of dimension d > 4 to the Lagrangian. 
These non-renormalizable operators are made out of the SM flelds, and invariant under the 
SM gauge group [l-3]. They parameterize the effects of the high-energy degrees of freedom 
on the low-energy theory order by order. In principle, the operator coefficients are weighted 
by inverse powers of the scale of new physics Anp: 



SM 



d=5 
cff 



nc>d=G I 



^cff 



With ^^^^ 



(1) 



The only possible dimension- five operator, namely, ^Sf^g ^, which violates lepton number by 
two units, is the famous Weinberg operator [l| 



(2) 



which leads, after electroweak symmetry breaking (EWSB), to Majorana masses for the 
neutrinos. Here L and stand for the Standard Model lepton doublets and the Higgs fleld, 
respectively. At tree level, O^r can only be mediated by a singlet fermion, a tri plet scalar. 



or a triplet fermion, leading to the famous ty pe I 



mechanism, respectively (see also Refs. [15|, [l6|). Compared to the electroweak scale, the 



type II 



13l |. or type III [ij] seesaw 



masses of the neutrinos in all three cases appear suppressed by a factor v/A^^p, where v/ \/2 
is the vacuum expectation value (VEV) of the Higgs fleld. Substituting typical values, one 
obtains that the original seesaw mechanisms point towards the GUT scale. 



^ Note that we will distinguish between the two concepts of non-standard four-fcrmion interactions (NSIs) 
and non-unitarity (NU). These two concepts will collectively be called non-standard effects. 
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Except for neutrino masses, the dimension- six operators potentially affecting neutrino 
oscillations are, for non-standard four-fermion interactions (NSIs), operators of the types 



0^={EE){LL), {LL){LL), 
and, for non-unitarity (NU) coming from heavy singlet fermions, of the type 

= {L<f) i^ (0tL) , 



(3) 



(4) 



where we omitted flavor, spin, and gauge indices. Some of these effective operators result in 
corrections to the low-energy SM parameters and in exotic couplings. For instance, Eq. (jl]) 



izmg 



implies a correction to the neutrino kinetic energy. After re-diagonalizing and re-norma 
the neutrino kinetic terms jl?! . 18|. a non-unitary leptonic mixing matrix appears 15|, Il9 |. 



On the other hand, Eq. Q typically leads to lepton-fiavor-violating processes. 

In ordinary seesaw models, the operators generating neutrino masses [such as Eq. ([2])] and 
non-standard effects [c/., Eqs. ([3]) and (jlj)] are both mediated by the same heavy seesaw par- 
ticles, and therefore, there might be connections between them in some cases, in particular, 
in the type-II seesaw. However, non-standard effects are usually suppressed dramatically 
by the scale of the seesaw threshold, which is typically not far away from the GUT scale. 
In some low-scale seesaw models, the smallness of neutrino masses is protected by other 
suppression mechanisms rather than the GUT scale, such as radiative generation 20|-|31|. 



32h49| . or neutrino masses from higher-than-dimension-five 



small lepton number breaking 

nn 

50H56[. In these tors may be introduced at the TeV 



effective operators 

scale, which typically lead to observable NU effects as well. As another example, in Table 2 



of Ref. 561], a number of possibilities to generate small neutrino masses together with NU 
are listed, where the neutrino mass originates from a dimension- seven operator. Of course, 
the heavy seesaw particles may also be directly searched for at future colliders, in particular, 
at the Large Hadron Collider via the lepton-flavor- violating decays 57f |. 

The NSI operators in Eq. ([3]) are typically connected to the charged lepton flavor violation 
by SU(2) gauge invariance, and constrained by lepton universality tests. This implies that 
stringent bounds exist for the NSIs (we list the current bounds in Sec. IIIIDI ). On the other 
hand, if one has a theory for which the (gauge invariant) NSIs in Eq. ([3]) appear without 
charged lepton flavor violation, the NSIs present in neutrino production at a neutrino factory 
and the neutrino propagation are correlated in a particular way, and the flavor structure 



of the NSIs is strongly constrained [58|]. Therefore, we focus on this class of operators in 
the following. The most prominent example for such a theory is the exchange of a heavy 
charged SU(2) singlet scalar, leading to an effective operator of the type 

(I^-L) (I-L^) , (5) 

where the dot denotes the SU(2) invariant product, i.e., ir^. This operator is antisymmetric 
in the flavor indices, and does not lead to charged lepton flavor violation. In fact, Eqs. (jl]) 
and ([5]) are the only two dimension-six operators which lead to non-standard effects without 
charged lepton flavor violation.^ In the remaining parts of this work, we will use NSIs and 
NU to denote the non-standard effects coming from Eqs. and (jl]), respectively. 

In fact, the above mentioned correlation between source and matter NSIs at the neutrino 
factory increases the experimental sensitivity to the NSI parameters dramatically 59|, sim- 
ilar to the NU case 60|. Furthermore, it is well known that NU can be re-parameterized in 
terms of NSIs 611] • We will demonstrate that for a neutrino factory, the mentioned NSIs are 
phenomenologically very similar to NU, but for completely different fundamental reasons, 
which makes it hard to distinguish them. On the other hand, at tree level, the natural im- 
plementation for the NSIs are the scalar bosons leading to the operator in Eq. ([5]), whereas 
the NU is mediated by SM singlet fermions. Therefore, distinguishing between the NSI and 
NU operators without charged lepton flavor violation is basically equivalent to differentiat- 
ing between scalars and fermions as mediators, at least to leading order at tree level, and 
therefore theoretically very interesting. Note that, besides leptonic NSIs, there may exist 
non-standard neutrino-quark interactions stemming from some Grand Unified or i?-parity 
violating supersymmetric theories which could also affect neutrino oscillations, but such 
scenarios lie beyond the scope of this work. 

The NSIs and NU have been extensively studied in the literature, from both theoretical 
and phenomenological point of view. In particular, it has been pointed out that the sub- 



leading effects generated by NSIs |62l-l65| add to the standard matter effect [66 



67] and also 



introduce new sources of CP violation; then a neutrino factory or a superbeam experiments 



^ In principle, there could be four-fcrmion interactions with two quarks and two Icptons contributing to 
neutrino detection or matter effects at a neutrino factory. However, it can be shown that these cannot 
produce NSIs without charged lepton flavor violation if written in a gauge-invariant way, because the 
up and down quark contributions carry different coefficients which makes it impc^sible to cancel both 
contributions simultaneously (without canceling the NSIs). See, for example, Ref. [l22| for details. 
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are adequate places to study their effects [68 
in the NU framework are quite similar, the same future 
(or measure) the additional rotations and phases 19 



. Since the sub-leading effects presented 



46 



aci 



47 



ities can 



De used to constrain 



60 



61 



76 



84|. 



In this work, we are mainly interested in studying the experimental signatures from 
both types of non-standard operators and trying to understand whether it is possible to 
disentangle them by using a neutrino factory facility or not. To this end, we will first present, 
in Sec. [Til the general formalism depicting neutrino oscillations in matter with both source 
and detector effects being included. These formulas are model-independent, and could be 
used in theories with both a unitary leptonic mixing matrix or a non-unitary one. In Sec. IIIH 
we will then classify the higher-dimensional operators according to the mediators, and figure 
out the NSI effects induced by different non-renormalizable operators. We also summarize 
the current bounds on the non-standard parameters discussed in this work. In Sec. IIV| we 
briefly discuss several possibilities how to determine the origin of the non-standard effects. 
Section |V] is devoted to a detailed analytical discussion of the transition probabilities useful 
for our analysis, as well as to the presentation of our simulation techniques and the numerical 
results to show the prospects of searching for the origin of NSIs in a future neutrino factory. 
Finally, a brief summary is given in Sec. I VII 



II. NEUTRINO OSCILLATIONS WITH NON-STANDARD EFFECTS 

In this section, we describe neutrino oscillations with non-standard effects from heavy 
mediators. Note that, NU essentially affects the couplings of neutrinos to gauge bosons W 
and Z, which are actually integrated out when describing neutrino oscillations, in particular, 
in the presence of matter effects. Effectively, NU leads to four-fermion interactions similar 
to these from NSIs, and therefore, both NSIs and NU can be described by using the same 
parametrization. For this part, it will be useful to treat both classes within the NSI frame- 
work, since the source, propagation, and detection effects are a priori treated independently 
for NSIs. Therefore, we use the NSI parametrization in this section. Namely, we present the 
oscillation probabilities including any type of non-standard effects in term of NSI parameters 
Eajs in despite of their possible origins. 

In order to perform more precision measurements on neutrino mixing parameters, an 
intense high-energy neutrino source together with a long-baseline setup is proved to be the 
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best choice 85|. Similar to the standard matter effects in long-basehne neutrino oscillation 
experiments, NSIs can affect the neutrino propagation by coherent forward scattering in 
Earth matter. In the language of effective Hamiltonians, the time evolution of neutrino 
flavor eigenstates in the presence of NSIs is described by 



NSI ; 



(6) 



where 



"Ho = -^Udi&g{m\,ml,ml)U'^ 
-Hm = diag(K:c,0,0), 
"Hnsi = Vcc^"^ , 



(7) 
(8) 
(9) 



with Vbc ^ V2GFNe arising from coherent forward scattering and denoting the electron 
number density along the neutrino trajectory in Earth. The vacuum leptonic mixing matrix 
U is usually parameterized in the standard form by using three mixing angles and one CP 
violating phase 



U 



C12C13 



S12C13 



"■S12C23 — Ci2S23Si3e"^ C12C23 — Si2S23Si3e"^ 



■523C13 



\ 



■5l2'523 ~ Ci2C23Si3e'^ — C12S23 — Si2C23Si3e'^ C23C13 



(10) 



where Cij = cos6ij, Sij = sinOij (for ij = 12, 13, 23), and 6 is the Dirac CP-violating phase. 



In analogy to the vacuum Hamiltonian "Ho in Eq. 
be diagonalized through a unitary transformation 



the effective Hamiltonian "H can also 



1 



n = ^^diag {mj 



~ 2 ~ 2 ~ 2 



(11) 



where m? 



1,2,3) denote the effective mass squared eigenvalues of neutrinos and U is 



the effective leptonic mixing matrix in matter. Note that, in writing down Eq. ( |TT|) . we have 
aheady taken into account the Hermitian property of "H. The explicit expressions for U and 
mf can be found in Ref. 



In addition to propagation in matter, production or detection processes can be affected 
by NSIs. The neutrino states produced in a source and observed at a detector can be treated 
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as superpositions of pure orthonormal flavor states 65|, |87 



K) = K) + = a + ^iu\ur^) > (12) 

l3=e,iJ,,T 

where the superscripts 's' and 'c?' denote the source and the detector, respectively. Note that 
the two flavor indices follow from the postulate of the coherent contribution to the source 
or detection effects, and more generally, there could be incoherent contributions, which will 



not be considered further 



87 



For NSIs, the parameters at sources and detectors are not necessarily correlated. Only 
if the production and the detection are exactly the same process with the same other par- 
ticipating fermions (e.g., beta decay and inverse beta decay), the same quantity enters as 



811 • l^ is important to keep in mind that these NSI parameters are experiment- 



and process-dependent quantities. In the following, we will mainly focus on the source and 
detector NSIs defined for a neutrino factory and specific processes. For NU, however, source, 
propagation, and detection effects are correlated in a particular way, as we will discuss in 
the next section. 

Including all the NSI effects into the neutrino oscillations, we arrive at the amplitude for 
the process ui 



a-y 



\T AT 



A' il + e' 



[A + e'A + Ae"^ + e'Ae'^] 



(14) 



where L is the propagation distance and A is a coherent sum over the contributions of all 
the mass eigenstates z/. 



Ar 



(15) 



With the above definitions, the oscillation probability is given by 90 1 

PK^4) = 1^-/3(^)1' 



+2 5^Im(j:,J-i;)sin^'"^^'^ 

i>j 



2 Am|L 
AE 



2E 



(16) 



where 



7 7 7,P 

A salient feature of Eq. f lT6|) is that, when a 7^ /3, the first term in Eq. f|T6|) is, in general, 
not vanishing, and therefore, a flavor transition would already happen at the source even 



before the oscillation process and is known as the zero-distance effect [9l|. Although the 
effective leptonic mixing matrix in matter U is still unitary, the presences of NSIs in the 
source and the detector prevent us from defining a unique CP invariant quantity like the 



standard Jarlskog invariant j92| . New CP non-conservation terms, which are proportional to 
the NSI parameters and have different dependences on the quantity L/ will appear in the 
oscillation probability. Another peculiar feature in the survival probability is that, in the 
case of a = /3, CP-violating terms in the last line of Eq. f ll6p may, in principle, not vanish. 
Note that Eq. f lTB]) is also valid in the case of a non-unitary leptonic mixing matrix. In 
the minimal unitarity violation scheme, the NU effects are parameterized by using similar 
e parameters as in the case of the source and detector NSI effects, but with the relation 

In what follows, we will discuss the possible origin of NSIs in Eq. ( !T6l) together with the 
correlations among the NSI parameters. 



III. ORIGIN OF NON-STANDARD EFFECTS 



As mentioned in the Introduction, non-standard effects naturally emerge from most fun- 
damental theories beyond the SM, and can, in general, be described by a series of higher- 
dimensional non-renormalizable operators after integrating out the heavy degrees of freedom 
in the underlying theory. The only dimension-five {O^) operator is the well-known Weinberg 
operator in Eq. ([2]), which gives birth to the masses of light neutrinos. For dimension-six 
operators (O^), depending on different mediators, there are, in general, two different kinds 
of non-renormalizable operators responsible for non-standard effects at tree-level. For scalar 
mediated dimension-six operators , four-fermion interactions are involved, which usually 
break the lepton flavor (and lepton universality), but conserve the unitarity of the leptonic 
mixing matrix. On the other hand, fermion mediated dimension-six operators O-^ correct 
the kinetic energy terms of light neutrinos, which violate the unitarity of the leptonic mixing 

8 



matrix as a consequence of the mixing between light and heavy neutral fermions. As for 
dimension-eight or higher operators, both two types of non-standard effects can be induced. 
In the following, we will discuss the possible non-standard effects stemming from different 
kinds of higher-dimensional operators. 



A. Dimension-six operators mediated by scalars 

If new scalars are introduced, the (leptonic) dimension-six NSI operators mediated by 



these at tree level below the EWSB scale are usually given by 



0^ = 2^/2Gf (e^/«)°J {v^YVlv^) {1'YVl,rI,) , (18) 

where £ denote the charged leptons. Here Gi? is the Fermi coupling constant and Pl and 
Pr are the left- and right-handed (chiral) projection operators, respectively. Note that, 
in writing down Eq. (fT8|) . we do not require gauge invariance. If SU(2) gauge invariance is 
imposed at the effective operator level and it is required that all the charged-lepton processes 
vanish, only the NSI operators of left-handed fields [the second type in Eq. ([3])] survive, 
which are antisymmetric in the flavor indices, i.e., a 7^ 7 and (3^5. Such operators can 
be naturally realized in theories with an SM SU(2) singlet singly charged scalar, leading to 



the operator in Eq. (|5]) 94l. |96|-|99|.'^ Therefore, the observation of a dimension-six operator 
of this type in the neutrino sector may be interpreted in terms of a heavy scalar boson 
as mediator. Note that, conversely, the consequence of any theory which does not lead to 
charged-lepton-fiavor-violation and produces dimension-six operators is, in general, that the 
antisymmetric conditions must hold 58|. 

For a neutrino factory, the (leptonic) NSI effects induced by Eq. ( !T8|) are relevant for the 
source, but not for the detector, since the detection processes involve quarks. Compared 



However, if appropriate cancellation conditions apply, one can also have theories without charged lepton 
flavor violation with triplet scalars, and singlet and triplet vectors. This can be read off from Table 2 in 



Ref . [58| : the effective operators generated by the exchange of the different mediators have to be combined 
such that C\]^ = — C£^. In addition, there may be loop- induced dimension-six operators. Since in these 
cases the theoretical interpretation is obscure (there maybe other effects induced by these operators), we 
focus on the singlet scalar interpretation in the following, which is the simplest one. 
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with Eqs. ( !T2|) and ( !T3|) . one can easily read ofF^ 



4" = i^X. , (19) 



'e/3 — 

-NF _ ( cL\e^^ 



<J = (^")«^ • (20) 



Here the dominating effect comes from the left-handed component at the production process 
due to the helicity suppression of the right-handed component. Note that we use the label 
"NF" to mark that these NSIs are (production) process dependent quantities only relevant 
for a neutrino factory, whereas potential NSIs at a superbeam are, in general, completely 
uncorrelated. In addition, note that by giving up two of the four flavor indices, these 
parameters violate CP and even CPT explicitly. For instance, e^^ may have some interesting 
effects, while e^J is completely irrelevant, since the beam does not contain any v^-. 
The leptonic NSI effects in matter are only sensitive to the vector component as 

ej. = {e'^n: + , (21) 

where two charged leptons are restricted to electrons. From Eqs. ( IT8|) and ( 12T]) . we find that 
e^p = {e^^y. In addition, since the beam typically transverses ordinary matter consisting 
of electrons, these parameters are (almost) experiment independent. 

Now, if we apply the antisymmetric condition (from gauge invariance and no charged- 
lepton-fiavor- violation), the matter NSI parameters (for a = e,fi,T) and the source NSI 
parameters e^^ and e^f are forbidden. In addition, we have the following relations {sSj 

^™ = -«")* , (23) 

with both and e^^ being uncorrelated. Of course, if charged lepton flavor violation is 
only suppressed, the relations in Eq. (l23l) only hold to some degree. However, we assume 
that the underlying theory does not produce charged-lepton-fiavor-violation, such as the 
mentioned singly charged scalar. 



^ Note that there is no standard ly^ production in a neutrino factory, and hence, there is no corresponding 
NSI parameter like e^^. 
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B. Dimension-six operators mediated by fermions 



In general, gauge invariant theories extended the SM with the tree-level exch ang e of heavy 
neutral fermions result in a dimension-six operator in the form of Eq. 1^ 19 1 



= c^p (lJ) i^ ( A ^ (24) 



with Ca/3 being the model dependent coefficients and being related to the Higgs doublet by 
= ir20*.^ After spontaneous breaking of the SM gauge symmetry, the operator defined in 
Eq. leads to a correction of the neutrino kinetic energy term, and hence, the leptonic 
mixing matrix deviates from unitarity. Note that the NU effects only make sense by means 
of effective theories, while unitarity will be restored once the "full" theory is taken into 
account. 

In the case of a non-unitary leptonic mixing matrix, the mass eigenstates (the phys- 
ical states) of neutrinos are linked to their flavor eigenstates by means of a non-unitary 



transformation 



84| 



\uf) = N\u^) = {l+ri)U\u„,), (25) 

where r] ^ —cv'^/2 is a Hermitian matrix and f/ is a unitary matrix diagonalizing the 
neutrino mass matrix. Note the similarity to Eqs. ( !T2|) and (fT3|) with respect to the source 
and detector effects; however, one can also see the difference compared to NSIs: for both 
NSIs and NU the matter effects are given in terms of lum), which, however, appear on the 
left-hand-side of Eq. f l2^ and implicitly in Eq. f[T^ . In the NU case, the flavor basis, through 
which the NC and CC interactions are deflned, is slightly shifted by rj. In the NSI case, 
additional contributions at source and detector may be present, which do not necessarily 
affect the properties of the weak interactions in matter (which are still deflned with respect 
to the original flavor eigenstates, and the link between mass and flavor eigenstates remains 
unitary). The time evolution of neutrino flavor eigenstates is given by 

i^l^/) = (1 + ^)^(1 + ^)"'!^/). (26) 



^ Apart from fcrmion singlets, this operator could also be realized in other frameworks once cancellations 
or loop- induced effects arc taken into account. 
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where 

H = Ho + (1 + v)^ diag (Vcc - Vnc, -Vnc, -Vnc) (1 + v) 

~ 7/0 + + {(-H™ + H^c) , r/} + . (27) 

Here = —diag (Vnc, ^nc, ^nc) and Vnc — V^GpNn, with being the number density 
of neutrons in Earth matter. Since NU effects are sub-leading order effects, one can safely 
neglect the terms proportional to 77^ in Eq. fl27|l . The pure NC contribution [the third term 
in the second row of Eq. (127|) ] is flavor blind, and hence can also be dropped. Then, by 



comparing Eq. f l27p with Eq. (E]), we obtain the relations 98 1 



■Hnsi = {{nm + nNc),v} , (28) 



^aB = Vae^fje + VeP^ae " ^^VaP , (29) 

for neutrino propagation in matter, and 

= = VafS , (30) 

for the source and detector effects. Note that, for neutrino propagation in realistic Earth 
matter, A'^e — holds to a very good precision. Therefore, we have approximately 

^ee — 2?7ee 5 — ~Vuii- ' ^/xt — ~V^lt 5 ^tt — ~Vtt , (31) 

ations only hold up to the percentage 



. This implies that and are 



together with ~ = 0. In practice, these cancel 
level, depending on the composition of the material 
expected to be suppressed by two orders of magnitude compared to the other parameters. 

In comparison with the NSI effects in the previous subsection, a salient feature is that, 
both the source and detector effects exist, and they are process and experiment independent. 
In the mean while, they always lead to interferences between the non-standard effects and 
the standard oscillation effects. 

C. Other categories of non-standard effects 

Apart from purely leptonic NSIs, there could be NSI operators involving quarks, or 
NSIs from leptonic dimension-six operators, which do lead to charged-lepton-fiavor- violation. 
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Furthermore, NSIs may be induced by dimension-eight or higher operators. For dimension- 
eight operators at tree level, however, effects (either NSIs or NU or both) are induced 



as we 
tests 



1, or exotic fermions appear, which are strongly constrained by electroweak precision 



58 



98 



There exists the principle possibility that the dimension-six NSI operators 



coming from different sources cancel and the leptonic NSI effects originate exclusively from 



dimension-eight or higher operators [58|. In such a case, the NSI operators might not 
produce charged lepton flavor violation either. The NSI operators involving four lepton 
doublets [second category in Eq. ([3])] allow then for a connection between source and matter 
effects, which could be different from the one discussed in Sec. IIII A[ On the other hand, 
if NSIs come from operators with two lepton doublets and two singlets [first category in 



Eq. Q], only matter NSIs will be induced 



In summary, there is a third category of non-standard effects, for which the source, 
detector, and matter interactions may be uncorrelated or correlated in a different way than 
in the previous subsections. However, these possibilities are either suppressed by higher 
orders of the new physics scale (dimension-eight operators), or face other constraints. For 
the sake of simplicity, we do not discuss these categories in detail any further, but we will in 
some cases point out when observations correspond to this category "Other". One should 
also keep in mind that data can only be interpreted in certain ways induced by new physics. 
We mainly discuss the interpretation in terms of dimension-six operators. 



D. Bounds on the dimension-six operators 

The bounds on the dimension-six operators are, in fact, model dependent. The experi- 
mental bounds mainly come from the lepton flavor violating decays £q, — )■ £^7, the univer- 
sality test of weak interactions and the invisible decay width of the Z-boson and have been 



summarized in Ref. 98||. For example, for a scalar mediated operator, one has 



k;;i = k^/l = d < 8.2xl0-\ (32) 
= le;:^ I < 1.9xl0-^ (33) 



'TT I 



< 8.4 X 10"^ , (34) 
le^/l < 7.5 X 10-2. (35) 
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For the non-unitarity effects induced by the O'^ operator, if the mediators are heavier than 
the electroweak scale, one has upper bounds on the 77 parameters^ 

^ 2.0 X 10^3 6.0 X 10~^ 1.6 X 10"^^ 
8.0 X 10-^ 1.1 X 10-3 
~ ~ 2.7 X 10^3 



\r]\ < 



\ 



(36) 



If the mediators are hghter than the electroweak scale but above a few GeV, the above 
bounds still apply except |?7e^| < 9.0 x 10^"^ has to be employed because of the restoration of 
unitarity in the Z-decay. In the case that NSIs come exclusively from d > 8 NSI operators, 
the severe constraints from universality test may not apply coherently, and hence, the bounds 
on NSI parameters are rather loose. 

The more gene ral, model independent NSI bounds for a neutrino factory (at 90 % C.L.) 



are given by 



100| 



^^0.025 0.030 0.030^ 



0.025 0.030 0.030 
0.025 0.030 0.030 



(37) 



and 



0.06 0.10 0.4 
0.03 0.10 
~ 0.16 



I mi? I ^ 



0.14 0.10 0.27 
~ 0.03 0.10 
~ ~ 0.4 



(38) 



for left- and right-handed NSIs, respectively. 

From this comparison, one can already see one dilemma which could be called the "NSI 
hierarchy problem" . While the model-independent bounds are relatively weak compared to 
the bounds on the dimension-six operators, the theory leading to such large non-standard 
effects cannot be that straightforward. For instance, if the NSIs came from dimension- 
eight operators, they would be naturally expected to be of the order w^/A^p ^ 10"^ for 
Anp = 0{1) TeV (as expected by the hierarchy problem). However, in this case, one cannot 
use particular correlations between source and matter NSIs to enhance the sensitivity. This 
means that the non-standard effects may in either case be on the edge of the sensitivity of 
a neutrino factory (see, e.g., Ref. [3] for matter NSIs). 



^ Note that, compared with the bounds on NN^ m Ref. ^9^, the constraints on ij are strengthened by a 
factor 2, since NN^ ~ 1 + 2r/ according to Eq. 
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TABLE I: Allowed parameters from the discussed dimension-six effective operator classes in a 
neutrino factory (i' -factory) and a superbeam experiment (SB). 

IV. TESTING THE ORIGIN OF NON-STANDARD EFFECTS 

In this section, we qualitatively discuss how the origin of the non-standard effect can be 
determined. In the rest of this work, we then focus on one particular example - at a 
neutrino factory. 

As we have shown, no matter of their possible origins, the non-standard effects can 
always be re-parametrized in terms of e'^, e'^, and e"^. Therefore, we have treated them 
as independent parameters in Sec. [ITl which can be used for any category. However, for 
a given experiment (such as a neutrino factory), e"^ will be a particular function of e'^ in 
different frameworks, and not all the NSI parameters are allowed [c/., discussion around 
Eqs. fl22|) . fl23|) . and fl3T]) ]. We summarize in Tab. [T] which non-standard effects are allowed 
for a neutrino factory and a superbeam if the origin are the discussed leptonic dimension-six 
operators. 

At this point, we want to emphasize again that the relationships between source and 
matter effects in Eqs. (!22|) . (!23l) . and (I3T1) are very similar, especially for the effects easiest 
to access. Consider, for instance, 

e'i^r = -(^f )* (NSIs) , (39) 

^;r. = -V,r = -4. (NU) . (40) 

The origin of these relationships is very different. The NSI relationship relies on the degree 
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that charged lepton flavor violation is suppressed, whereas the NU relationship relies on the 
equality of the CC and NC potentials in Eq. fl29|l . 

If one wants to distinguish the origin of the non-standard effects, it is not sufficient to 
determine the phase of e^^, since one can always fit the data with NSls with one phase or NU 
with minus this phase. In principle, one needs an independent test of and 5™ (including 
phases). However, we will demonstrate that the equality between source and detector effects 
for NU leads to additional differences between the effects. 

From the previous discussions, we have qualitatively four possibilities to determine the 
origin of the non-standard effects (c/.. Tab. [T]): 

• Distinguish by appearance of certain parameters. If, for instance, is found, it 
cannot be interpreted as O'^ , but as O"^ (or the category "Other"). If, on the other 
hand, is found, it has to come from the category "Other", such as a higher- 



dimensional operator. See, e.g., Ref. [75|, for the bounds expected for these parameters 
at a neutrino factory. While is one of the most discussed NSI parameters in the 
literature and can be relatively well constrained, adds to the standard matter 
effect and the bounds are limited to the precision by which the matter density profile 



is known 



75| . Another such parameter is e™, which can only come from "Other", and 



is discussed in Ref. [73|. In addition, e^^ or could discriminate between and 
O'^. At a neutrino factory near detector, this measurement will be limited by charge 
identification. 

Distinguish by bounds (c/.. Sec. IIII Dp . If large enough effects beyond the dimension- 
six operator bounds are found, but below the generic bounds, they have to come from 
the category "Other". If, for instance, e^^ ~ 10"^ is measured, it could come from 
O^, but is excluded for O^. 

Distinguish by experiment class. Since the source NSIs are production process- 
dependent parameters and the NU parameters are fundamental, using a different ex- 
periment class can help to disentangle the effects. For instance, if e^^ is found at the 
near detector of a neutrino factory, it may come from or O^. If it is a fundamental 
parameter from O-^, it has to appear at a corresponding superbeam near detector, 
such as the anticipated Ma in In jector Non-Standard Interactions Search (MINSIS) 



project at Fermilab (USA) jl01| . Otherwise, C-^ could be excluded. Therefore, in 
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order to disentangle pro cess- dependent NSIs from fundamental NU, neutrino factory 
and superbeam near detectors are complementary. 

• Distinguish by correlations. For example, = = — holds for O-^, while 
e^^ = —e"^* together with = holds for . The other correlations in Eqs. f l22|) 
and ( 13T]) will be less accessible at a neutrino factory, since either there are no Uj. in the 
beam or the measurement of el^ and e*^ will be intimately connected to the knowledge 
of cross sections and fluxes. 

Since especially the last category requires some more detailed understanding, we focus on e^r 
for the rest of this study. In this case, it has been pointed out in the literature that the U/i ~ ^ 
Ut- a nd Ufj^ — z/^ channels provide us with the best sensitivities (see, e.g., Refs. 



61 



82 



84 



102 



1061], where the phenomenological importance of the fifi and fir transitions for searching new 
physics has been stressed. Therefore, we will mostly concentrate on these two channel in 
the subsequent analysis. 



V. NON-STANDARD EFFECTS IN THE /i-r-SECTOR AT A NEUTRINO FAC- 
TORY 

Here we focus on the measurement of e^r, the differences between the NSI and NU 
sensitivities, and the ability to determine the origin of the non-standard effects at a neutrino 
factory. 



A. Analytical considerations 

In this section, we discuss some of the relevant analytical properties of the transition 
probabilities useful to obtain knowledge on the output of our numerical simulations. We 
will concentrate on the appearance probability -P(t'^ — z^^) as well as the survival probability 
P(z/^ — )■ u"^), keeping all ^^'J"'"^ but e^'™''^ vanishing. Note that, since we are only considering 
neutrino factory setups, we will drop the upper NF label on the non-standard parameters. 
The expressions for these probabilities can be obtained by applying the general formula 
given in Eq. f|T6l) : to further simplify the results, we consistently expand up to first order in 
the small quantities 6*13, Aml-^/ Aml-^^, and e^/j. Keeping all the source, matter, and detector 
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effects, the survival probability reads 



P, 



MM 



sin 26*23 sin 



- (Ree^^ + Ree^^) sin 4^23 sin^ 



AL 



+ {imef^ + Ime^^) sin 26*23 sin 

aL 2 

sin 26^23 sin 



-Re sin 2^23 



2E 



AL 
'2E 
/AL 
V 2E 



4a 
'A 



cos 26*23 sin 



AL 

Ie 



(41) 



where we defined a = 2EVqc and A = Am|^. This result agrees with Eq. (35) in Ref. 81 |. 
Note that the presence of an explicit CP-violating term in the disappearance channel is 
not a signal of breaking the CPT symmetry, but the result of considering different initial 
and final states. This "CPT violation" is eliminated once the considered production and 
detection processes are the same with the same participating fermion, for which the relation 



Im£ 



Ime*^ holds. It is now straightforward to obtain the z^^i — t- transitions in the 



case of scalar and fermion mediated dimension-six operators. In the first case, we need to 
drop the dependence on in Eq. ( HT]) and use the mapping in Eq. ( 12^ : choosing the source 
parameter as the relevant parameter, we obtain 

'AL^ 



F. 



1 — sin 26*23 sin 



4E 



—Re sin 46*23 sin 



AL 

Ie 



+Re£^^ sin 2^23 



aL 
2E 



sin 26*23 sin 



Imef,^ sin 26*23 sin 
4a 



flT 

AL 



AL 



2E 



A 



2E 

cos^ 26*23 sin^ 



m 



(42) 



In the case of fermion mediated operators O^, we can use Eqs. f l30|) and fl3T]) to obtain 

/AL' 



P. 



T 



sin2 2^23 sin^ 



\4E 



-2 Re sin 46*23 sin^ 



(AL 
V4E 



+Ree" sin 2B 



23 



aL 2 

— sm 26*23 sm 
2E 



AL 
2E 



4a 2 /I 2 
— cos 26*23 sm 
A 



/AL 

[Je 



(43) 



Equations. fH2|) and f H3|) have exactly the same form, except for the second lines, which are 
different from each other. In Eq. ( H3l) . compared to Eq. ( |4T]) . it is clear that a cancellation 
between the imaginary parts of the source and detector parameters as well as a sum of their 
real parts is at work. This is not an accidental fact, but can be proved to happen for any 
disappearance probability; in fact, from the general expression in Eq. ( IT6|) . the true CP 
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violating terms appearing in the last line of Eq. (fT6|) disappear once we use the definition 
for J^p in Eq. f lT7|) computed for the same flavor a = (3 and apply the mapping relations 
given in Eqs. (15U]) and CT. 

For the numerical simulations to follow, it is useful to adopt a short baseline expansion 
for the z/^ — t- Vj. transition, namely up to second order in small e's and flrst order in the 
quantity AL/AE. For the general NSI effects, the transition probability is given by: 

P,. = sin22^23 (|^)' + l4.r + l<.r-2 (Imej^ + lme^,) sin2^23 

+2 Re (ej^ep) (44) 

and it reduces to the following expressions for the scalar and fermion mediated operators, 
respectively: 

P;f, = sin2 2^23 (II) +|<j2-2Im5^^sin2^23 (||) , (45) 

P^r = «in^2^,3 (||)' + 4|<J^-4Im.^,sin2^23 (^) • (46) 

It is noteworthy that, in these two equations, the third term dominates over the second if 
< AL/{2E) ~ 10^^ at the energy threshold (about 1 GeV) of a neutrino factory for 
L ~ 1km, whereas the flrst term is suppressed by L^. This means that the non-standard 
CP violation from Ime^^ is, in principle, measurable in a near detector for large enough 
statistics. The typical (relative) statistical error for the near detectors considered is about 
10~^ to 10"^, which is to be compared with {e^^^)"^. In conclusion, the currently considered 
detectors are on the edge of measuring that effect. 



B. Simulation techniques 



We continue to the numerical simulations of searching for different NSI effects at a 
neutrino factory. In our analysis, we mostly follow the International Design Study (IDS- 
NF) 1071 baseline setup, which consists of 2.5 x 10^° useful muon decays per polarity and 
year with the parent muon energy i?^ = 25 GeV. The total running time is assumed to 
be ten years. Two magnetized iron calorimeters are assumed at L = 4000 km (flducial 
mass 100 kton) and L = 7500 km (fldu cial mass 50 kton), respectively. The description 



of the neutrino factory is based on Refs. 



107Hll0j . In addition, we consider in some cases 
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OPERA-inspired (magnetized) Emulsion Cl oud Cham ber (ECC) near i/^-detectors. The 



z/g — )■ channel description is based on Refs. 
have the same characteristics as in Ref. 



109 



. The — )■ v- r cha nnel is assumed to 



governed by Refs. [6ll. Il02| . Since we assume 



that the hadronic decay channels of the r can be used as well, we assume a factor of five 
higher signal and background than in Ref. i.e., 48 % detection efficiency. Because we 

are mostly interested in the — )■ channel (for which the Ve channel is only a small 
perturbation in the presence of only), we add the and Vr events as a conservative 
estimate (however, there is little impact of this assumption for the parameters considered). 



The near detector treatment is based on Ref. [59| with respect to geometric effects of de- 
cay straight and detector geometry. For the high-energy neutrino factory, we consider the 
following near detectors, based on the simulation described above: 

ND-L Large (OPERA-like) size, fiducial mass 2 kt, c? = 1 km (distance to end of decay 
straight), 

ND-M Medium size (e.g'., SciBar-sized), fiducial mass 25 t, = 80 m, 
ND-S Small size (e.g'., silicon vertex-sized), fiducial mass 100 kg, = 80 m. 



ON D@ 1 30km OPERA-like at intermediate baseline L = 130 km, as proposed in Refs. 61 



105| . in order to improve the sensitivity to the non-standard effects. 



Note that is yet unclear if an ECC can be operated as close as 1 km to a neutrino factory 
because of the high scanning load. Therefore, alternative technologies may be preferable, 
such as a silicon vertex detector. In this case, other challenges have to be approached, such 
as the background from anti-neutrino charm production. These issues are currently under 
discussion within the IDS-NF. 



In addition, to a high-energy neutrino 
neutrino factory based upon Refs. 
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actory, we consider a low-energy version of the 
113j. It has = 4.5 GeV and a magnetized Totally 
Active Scintillator Detector (TASD) with 20 kt fiducial mass at L = 1300 km. The total 
running time is assumed to be ten years with 7 x 10^" useful muon decays per polarity 
and year. This luminosity is higher by a factor of 2.8 than the one of the standard neutrino 
factory, since all muons are put in one storage ring (factor two) and the frontend is optimized 
in a different way, which leads to another 40 % increase of the luminosity 113| . Note that 
we include two types of backgrounds for the appearance channels, one which scales with 
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the disappearance rates (such as from charge mis-identification), and one which scales with 
the un-oscillated spectrum (such as from neu tral c urrent events), both at the level of 10~^ 
(which is a factor of two hi gher than in Ref. llSj). We also include the — )■ z/g channel, 



based upon Refs. 



109 



114l |. As near detectors, we consider the above mentioned OPERA- 



like detector (ND-L), and the small (ND-S) with the same characteristics at a distance of 
20 m fro m th e end of the decay straight. For the decay straight, a length of s = 200 m is 



assumed 



115l | (needed to compute the effective baseline as desc ribed in R ef. [59|) 



For the experiment simulation, we use the GLoBES software 
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systematics. For the oscillation parameters, we use (see, e.g., Refs. 



117| with user-defined 
119^ sin2 2^i3 = 0, 



1181 



sin^^ia = 0.3, sin^ ^23 = 0.5, Am^^ = 8.0 x 10^^ eV^, Am^^ = 2.5 x lO-^eV^, and a normal 



mass hierarchy, unless specified otherwise. We impose exte rnal error s on Am2x and 612 of 



3 % each, and we include a 2 % matter density uncertainty 



120, 



121 |. 



C. Numerical results 



In Fig. [H we illustrate the discovery reach of the IDS-NF neutrino factory to the parameter 
e*^ in the two scenarios of fermion (left panel) and scalar mediated (right panel) operators. 
We define the discovery reach as the values of (true) parameters and [where = 
le^^l exp(i0^^)], for which je^^l = can be excluded at 90 % C.L. The standard oscillation 
parameters are thereby marginalized over. As we discussed above, in the case of fermionic 
operators O"^, the cancellation between the imaginary parts in the survival probability leads 
to the far detectors being sensitive to the real part of only, as it can be clearly seen 
from Eq. (H2l) . Therefore, there is hardly sensitivity to the corresponding e parameters in 
the case of the CP- violating phase = ±7r/2, as shown by the solid (red) curve in the 
left plot of Fig. [U In order to increase the sensitivity to the CP- violating phase 0^^, it is 
useful to include the effects of near detectors capable for r identification. In fact, as we can 
see from Eq. P6l) . an explicit dependence on the imaginary part of appears in and, 
depending on the type of the near detector, this term could be relevant. Note also that the 
le^^P term helps in increasing the sensitivity to the absolute value of e^^. In this plot, we 
also compare three different scenarios, in which we combine the IDS-NF neutrino factory 
setup with the near detectors ND-S, ND-M, and ND-L. The best combination is using the 
ND-L detector, mainly due to the larger mass. As a rough estimate, a sensitivity 3 x 10~^ to 
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FIG. 1: The discovery reach of the IDS-NF neutrino factory to e^^^ originating from O (left) or 
O'^ (right) operators (on right-hand-side of curves, 90 % C.L.). The shadowed areas indicate the 
current experimental constraints on the corresponding dimension-six operators. In both plot the 
effects of adding different near detectors are also shown. 



\b^^^\ can be expected in the presence of a near i/.r-detector. Note also that a short distance 
z/T--detector {L ^ 100 km) does not help much in improving the sensitivity, and hence, we 
did not include the relative curve in the plot. For a scalar mediated NSI operator, the 
situation is a bit different, due to the no n- vanishing contribution from the imaginary part of 
e*.^, see Eq. ( I45l) . Although this term is dominated by the matter-induced one, its effect is 
already visible in the right plot of Fig. [H where the CP violating term is responsible for the 
asymmetric behavior of the sensitivity curves with respect to 0^.^ = 0, and better sensitivity 
for ~ ±90°. Similar to the case, a better sensitivity could be expected once the 
larger near detector is taken into account (the scenario label NF + ND-L in the plot). 

Compared to standard oscillation physics, where the statistics in the far detectors limit 
the performance, the size of the near detector is very important for non-standard effects. A 
meaningful question could be how well the non-standard effects can be measured if they are 
not vanishing. An example is given in Fig. [21 in which the best-fit contours for the chosen 
NSI parameter = 0.001 exp(i7r/4) are plotted for both type of operators in three different 
situations, where the IDS-NF neutrino factory setup is alone or accompanied by one or two 
OPERA-like detectors at different baselines, namely at L = 1 km and L = 130 km. Using 
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FIG. 2: The 90 % C.L. contour plot of the N SI parameter e'l^ originating from O'^ (left plot) or 



(right plot), assuming the true paramet 



er e, 



0.001 exp(i7r/4). The analyzed experimental 



setups consider the IDS-NF neutrino factory alone and in combination with an OPERA-like near 
detector at a baseline of 1 km and/or 130 km. 



the other two options ND-S and ND-M do not improve the performance. Due to the fact 
that the far muon-detector is not sensitive to the imaginary part of (expect from a small 
effect for O'^), we observe in both panels that a standard neutrino factory without near 
detectors has almost no sensitivity to the CP-violating part sin0^^. The situation improves 
a lot in the presence of near detectors, especially when a short-distance detector located at 
L ~ 100 km is taken into account, in which the sin A terms in Eqs. f H5|) and fH6|) contribute 
to the appearance probability. In particular, the CP-violating phase of e*^ can be better 
reconstructed for the operator (left panel) because of the additional improving factor of 
two in the imaginary part term. 

Since the phenomenological signatures of O"^ and in terms of discovery potential and 
parameter measurements are quite similar, it is an important question to see if a neutrino 
factory-based setup is able to discriminate from O'^ so as to find hints on the origin 
of non-standard effects. We answer this question generating events by using the "true" 
parameters in the case of O"^ (C^*^)? and then fit the data with only (C^'^)- The results 
of such an analysis are the exclusion regions (right-hand side of the curves) shown in Fig. |3l 
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FIG. 3: Regions in the f]e^^|-(/)^^j-p/ane where the simulated e^^^ induced by one type of operator 
can be uniquely established, i.e., the other type of operator is excluded at the 90 % C.L. (regions 
on the right-hand side of the curves). Left panel: the simulated is induced by O-^ and fitted 
with . Right panel: the simulated e^^ is induced by and fitted with . In both panels, the 
discovery reach is also displayed. The experimental setup is the same as Fig. \^ 

where the dimension-six operators can be disentangled. For the IDS-NF neutrino factory 
combined with several different near detectors, the curves in the left panel show that there 
is just a very small region beyond the bound at the 90 % C.L. on the O'^ operators, where 
the data generated with (9-^ can be distinguished from the even if the OPERA-like near 
detector at the longer baseline is used. If is simulated, however, it can be distinguished 
from for a part of the parameter space beyond the current bound with ND-L, as shown in 
the right panel of Fig. [3], especially around the CP-conserving value = 0, ±7r, where the 
current experimental constraints on are not as stringent as for O-^ . The combinations 
with other near detectors are illustrated by the ND-S curves, where the sensitivity does not 
go beyond the current bounds. 

From Fig. [3l we can read off that there are substantial regions of the parameter space 
between the curves, where and can be distinguished at the neutrino factory itself, 
and the discovery reach of the neutrino factory, for which the neutrino factory will find a 
non-standard effect, but cannot classify it. These regions go beyond the current bounds. 
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FIG. 4: Same plots as Fig. [7] (left) and\^ (right), but for a low-energy neutrino factory alone and 
in combination with the ND-S and ND-S+ND-L. 

As mentioned before, an interesting discriminator might in this case be a superbeam v^-- 
detector, such as the MINSIS project. Such a detector could tell the fundamental effect 
from the process-dependent effect if the sensitivity is comparable to that of the neutrino 
factory, i.e., |e^^| ^ 10~^. For instance, if no effect is seen at MINSIS, but some effect is 
detected at the neutrino factory, it could come from , but not from O'^ . 

Recently, a low-energy neutrino factory [E^ = 4.5 GeV) has been attracting some atten- 
tion, for which the same kind of analysis in Figs. [1] and [3] can be repeated. In the left panel of 
Fig. m we show the discovery reach for the O'^ induced NSIs using the low-energy neutrino 
factory alone and in combination with different combinations of near detectors. We observe 
that, because of the r production threshold, there is no hope to search for NSIs originated 
from O'^ operators below the current experimental limits. Also, in the right plot of Fig. H] 
the same combination of experimental facilities is not able to exclude the operators at 
90 % C.L., since the exclusion regions are excluded by current limits already. Finally, in 
Fig. [5l we plot the CP discovery potential for both (left panel) and (right panel) 
induced CP violations. This is defined as the ensemble of true values of 0^^, which cannot 
be fitted with the CP-conserving values = 0, ivr at 90 % C.L. The combination of the 
standard IDS-NF neutrino factory with different large enough near detectors may discover 
CP violation, somewhat beyond the current bounds, especially for 0f,^ ~ ±7r/2. There are 
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FIG. 5: The 90 % C.L. of CP discovery potentials in the NU framework (left) and N SI framework 
(right). The facility combination is the same as Fig. O 



no qualitative differences between (9-^ and . 
VI. SUMMARY 

In this work, we have clarified the relationship between NSI and NU effects, where we 
have focused on non-standard effects coming from dimension-six effective operators when the 
heavy fields are integrated out. At tree level and without cancellations, these operators are 
mediated by scalars (NSIs) or fermions (NU), which means that the discrimination between 
NSIs and NU is interesting from a theoretical point of view, since it may reveal the nature 
of the heavy mediator. 

From the phenomenological point of view, the assumption of NSIs or NU, together with 
the assumptions of gauge invariance and vanishing charged lepton fiavor violation, has lead 
to particular correlations between source, propagation, and detector non-standard effects. 
These correlations can, in some cases, be used to disentangle NSIs from NU, relying on the 
measurement of individual parameters (such as for e™). However, for a neutrino factory, 
NSIs and NU look very similar for some parameters - for entirely different fundamental 
reasons. 

The most interesting case may be e^,-, which is, in principle, easy to find at the near de- 
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tector of a neutrino factory (there are no z/^ in the beam). However, the correlation between 
source and matter NSIs is basically the same for NSIs and NU. There is some discrimination 
potential coming from the fact that NU has modified source and detector effects, which, 
however, hardly exceeds the current bounds. Thus, the easiest way to discriminate NSIs 
from NU is the comparison with an experiment using a different neutrino production mech- 
anism, such as a superbeam. A possible project in that direction is the MINSIS project. In 
order to provide complementary information, a similar sensitivity is required, which should 
be significantly below 10^'^ for 

We conclude that differentiating between NSIs and NU should be one of the key priorities 
of searches for new physics effects, since the nature of the non-standard effect points towards 
the nature of the heavy mediator. The components necessary for this search are detection 
at least in near detectors, both at high-intensity superbeams and a neutrino factory. For the 
neutrino factory, a high enough muon energy is mandatory for the discussed non-standard 
effects searches, which means that the high-energy neutrino factory should at least be an 
upgrade option even for large 6'i3. In addition, for non-standard effect searches, the size 
of the near detector is very important, which means that for all applications, large enough 
detectors are needed. 
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